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Abstract 

This work explores a simple model of a slender, flexible structure in a uniform flow, 
providing analytical solutions for the translating, axially flowing equilibria of strings sub¬ 
jected to a uniform body force and drag forces linear in the velocities. The classical 
catenaries are extended to a five-parameter family of curves. A sixth parameter affects 
the tension in the curves. Generic configurations are planar, represented by a single first 
order equation for the tangential angle. The effects of varying parameters on represen¬ 
tative shapes, orbits in angle-curvature space, and stress distributions are shown. As 
limiting cases, the solutions include configurations corresponding to “lariat chains” and 
the towing, reeling, and sedimentation of flexible cables in a highly viscous fluid. Regions 
of parameter space corresponding to infinitely long, semi-infinite, and finite length curves 
are delineated. Almost all curves subtend an angle less than vr radians, but curious spe¬ 
cial cases with doubled or infinite range occur on the borders between regions. Separate 
transitions in the tension behavior, and counterintuitive results regarding finite towing 
tensions for infinitely long cables, are presented. Several physically inspired boundary 
value problems are solved and discussed. 


1 Introduction 

The catenary, or hanging chain, is one of the oldest problems in analytical mechanics. In the 
late 17th century, Hooke identihed its shape with that of an inverted arch. In keeping with 
scientihc customs of his time, he hid this knowledge in a very long anagram, appended as hller 
to an unrelated work on helioscopes [1, 2]. Not long after, in response to a challenge by Jacob 


*bchakrab@eng.ucsd.edu 
ihannaj@vt.edu 


1 



Bernoulli, the solution was found independently by Johann Bernoulli, Huygens, and Leibniz [3]. 
In the mid-19th century, the catenaries were recognized to be solutions of a broader dynamical 
problem, namely that of determining the shape equilibria of a moving string subjected to a 
uniform body force. This was perhaps hrst recognized by a Tripos examiner in 1854, one of 
whose victims, Routh, incorporated the fact in a mechanics text a few years later [4], By 
that time, the discovery had also been made by Airy (Astronomer Royal) and Thomson (Lord 
Kelvin), whose impetus to study this problem was the massive failure of the hrst attempt to 
lay transatlantic telegraph cable [5, 6]. Their treatments of the effects of the huid medium fell 
short, with Airy modeling drag as isotropic with respect to the cable geometry, and Kelvin 
neglecting it entirely. 

As shapes, the catenaries are a family of curves with a single scaling parameter that can 
effectively describe the length of a string hung, or arching, between two points. Additional axial 
how along the tangents of the curve serves only to increase the tension in the string, leaving 
shape equilibria intact. The tension increase has consequences for the stability of equilibria, as 
has been observed in translating belts [7] and in the “chain fountain” phenomenon [8, 9]. To 
the symmetry of the system under axial how, we may add another under uniform translations 
of the curve, which correspond to Galilean boosts that ahect neither the shape nor the tension. 

In this paper, we extend these classical results to include the ehects of linear, anisotropic 
drag forces, such as would arise if the translating, axially howing catenary were immersed in 
a highly viscous huid. The presence of huid breaks multiple symmetries and turns the one- 
parameter family of catenary curves into a hve-parameter family. The additional parameters 
beyond scaling are the magnitude and direction of the translational velocity, the magnitude of 
the axial velocity, and the anisotropy of the drag forces. We hnd that, in contrast to symmetry 
breaking by gyroscopic means [10, 11], some symmetry still remains, in that inertial forces 
continue to only ahect the string’s tension, not its shape. This is because Coriolis terms are 
absent in our problem; centripetal contributions to inertia behave identically in both cases. 

Slender bodies immersed in huid may be encountered across a wide range of scales in 
the industrial and natural worlds. High-altitude power devices and ohshore structures rely on 
pendant tethers, mooring lines, risers, and umbilical cables loaded by currents of air and water. 
Towed cable systems are widely used for retrieval and sensing under the ocean. Some aquatic 
and terrestrial plants bend and twist to avoid damage or uprooting, and the hexibility of such 
structures suggests adaptive motifs for drag reduction or control [12, 13]. Biohlm streamers 
form at high [14] and low [15, 16] Reynolds number. Chemical and huidic applications may 
employ microscopic wall-mounted tubes or rods in shear or other non-uniform hows [17-19]. 
Suspensions of hbers at low Reynolds number arise in classical and industrial sedimentation 
problems [20, 21]. Other examples include spinning yarn [22], bacterial hagella [23], and the 
ballooning threads of migrating spiders [24] . These systems often involve nonlinear drag forces, 
nonlocal ehects, nonuniform hows, and bending stihness, which are complications we do not 
address in this paper. Before tackling such problems, it behoves us to consider our simple 
system, and try to use our conclusions as a basis for generalization. 

Our augmented catenaries are generically planar curves, and admit an analytical descrip¬ 
tion as a hrst-order dynamical system for the tangential angle, with arc length as the in¬ 
dependent variable. This suggests an approach to organizing our results, by classifying the 
corresponding phase space in terms of the presence and location of hxed points and poles. We 
illustrate the role of the several parameters on the solutions. We also discuss the nature of the 
string tension, which, surprisingly, undergoes its own bifurcations in a manner noncoincident 
with those of the phase portraits. Our discovery and exposition of the solutions is aided by 
our initial disregard for any specihc physical boundary conditions or integral constraints, which 
are connected to the axes of our parameter space only through implicit integral relationships. 
However, towards the end of the paper, we discuss several boundary value problems in which 
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practical parameters like the positions of end points and the total length of the curve are 
specihed. 

The general solution we consider covers several interesting problems as limiting cases. 
Aside from the classical catenary, we may use our results to examine the towing of a neutrally 
buoyant cable (or, equivalently, such a cable tethered in a uniform flow), the reeling in and out of 
a belt or “lariat chain” [25] between two points, or the free sedimentation of a perfectly flexible 
hlament. One limit, that of a string snaking along its own tangents under the action of drag 
forces alone, is a singular axis of our parameter space, admitting only the very restricted set of 
straight-line solutions. In contrast, the hypersurface in parameter space where drag anisotropy 
vanishes is a degenerate case, admitting arbitrary shapes as freely sedimenting solutions. 

In Section 2, we pose our problem and provide the general form of its solutions. In Section 
3, we present examples of generic solutions and organize them according to their singularity 
structure in phase space. In Section 4, we present a few solutions of physically realistic boundary 
value problems. The literature on several related problems is briefly discussed in Section 5. The 
Appendices include brief details on the symmetry properties of the phase space, and an alternate 
derivation and representation of a subset of our solutions that includes the classical catenaries. 

Throughout, we will use the following terms interchangeably: “tension” and “stress”, 
“string” and “cable”. 

2 Formulation and general solution 

We hrst consider catenary curves in three dimensions in Section 2.1, before casting and solving 
the problem in two-dimensional form in Section 2.2. Additional three-dimensional aspects of 
the solution are reserved for Section 4.1. 

2.1 The three-dimensional problem 

Consider a curve X(s, t) parameterized by arc length s and time t, carrying an adapted or¬ 
thonormal frame ^t(s, t), h(s, t), b(s, whose spatial rotations are described by two extrinsic 
curvatures and T{s,t). Hence, 
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We seek solutions in the form of shapes translating with uniform velocity v, with a uniform 
axial flow of material given by a tangential velocity TcI^X = Tt. These velocities lead to linear, 
anisotropic drag forces with coefficients for drag along the tangent and i>j_ for drag in the 
directions of either of the two normals n or b. A uniform body force q acts on the string. 
These elements are shown schematically in Figure 1. 

Material time derivatives, taken at hxed s, provide the velocity and acceleration for the 
assumed form of motion. 


aiX = v + Tt, (2) 

a^X = T^dsi = T\n , (3) 

allowing us to write a force balance on the string, 

mT^dgi -|- i^ii (T -|- V ■ t) t -|- i/j_v ■ ^hh -|- bb j = dg (at) -|- q. (4) 
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The first term represents inertial forces (per nnit volnme) arising from centripetal accelerations 
of material moving along the rigidly translating shape. The mass coefficient m represents both 
the material of the string as well as added mass effects from the snrronnding flnid, while the 
body force q conld be linked to the density of the body, an electric held acting on a charged 
body, or any other nniform effects. The drag coefficients i/y and are taken to be constant 
over the range of velocities experienced by the body^. Finally, a is the tension in the string, 
which acts as a mnltiplier enforcing the inextensibility constraint = 1. 



Figure 1: A catenary X is an inextensible curve acted on by a uniform body force q. The 
body carries the orthonormal frame (t,n, b), and is translating with uniform velocity v, with 
a superimposed uniform material how Tt along its tangents. 

Equation (4) has several limiting cases of interest. When the drag terms vanish, we recover 
the classical catenaries. If, in addition, the body force vanishes, any arbitrary conhguration X 
with uniform tension a = mT'^ is a solution; these are the Routh “lariats” that snake along 
their own tangents [4, 29]. Arbitrary conhgurations X are also solutions when drag forces are 
present, but isotropic (z/y = z/j_ = z/); these sediment as rigid, tensionless bodies (T = a = 0) 
at a velocity such that z^v = q. 

As is true for many string problems, generic solutions of the equation of motion (4) are 
planar, meaning that r = 0 and b is a constant vector. This may be seen by considering the 
projections onto the frame vectors. 


z/|| (T -b V ■ t) = ds(T + q ■ t, 

(5) 

+ z/xv • n = (TK -b q ■ n, 

(6) 

z/xv • b = q • b , 

(7) 


taking repeated s-derivatives of (7), remembering that v and q are constant vectors, and using 
the relations (1). One hnds that the only fully three-dimensional solutions are the limiting cases 
already mentioned. Additionally, equations (5-6) may be used to show that the only lariats 
(v = q = 0, T 7^ 0) in the presence of tangential drag forces are straight lines (zt = 0). 

^In the Stokes flow limit, where the Reynolds number vanishes, the drag coefficients can be derived analyti¬ 
cally from the leading-order singular terms in a boundary-integral representation of the fluid velocity. They are 
such that 2z/|| = and are logarithmically singular in the filament radius (normalized by some larger cutoff 
length). The derivation is similar to those found in [26-28], and breaks down if the curvatures get too large. 
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Before discussing generic planar solutions, let us pause and consider these simplest of 
equilibria, straight lines, known to arise when a cable has one free end. For straight lines, the 
K terms in (6) are zero, and the distinction between any vectors orthogonal to the tangent 
disappears. An equation like (7) holds for any such vector, which implies that the straight 
cable lives in the plane spanned by any choice of nonvanishing, noncollinear velocity v and 
body force q. Imagine towing a cable with one end on a reel and the other end free. Given 
V and q, the angle of attack of the cable is determined, and (5) relates the tension a and the 
payout speed T. The tension is a linear function of s. If one end is free, the constant obtained 
by integrating (5) is used to set the tension to zero there. The payout speed does not affect the 
angle of attack, but is linearly related to the tension at the tow point, and there is a particular 
speed at which an untensioned cable (cr = 0) will be pulled off of the reel. For collinear v and 
q, parallel straight lines can move at any velocity, while for perpendicular straight lines there 
is a single velocity satisfying = q. Either case admits any value of axial flow T, although 
in the parallel case the distinction between v and Tt is ambiguous. Similar calculations can be 
performed for nonlinear drag forces, when the geometry is linear like this. 

2.2 The planar problem 

We may use the planarity of generic solutions to simplify their description. Note that while 
the curves X lie in a plane, this plane does not, in general, contain the body force or the 
direction of translation. Specihcation of the plane requires use of the binormal force balance 
(7) and an additional condition that will typically be provided by a boundary condition. We 
reserve discussion of this process until Section 4.1, where we discuss a boundary value problem 
for towing a cable with specified endpoints. For now, we note that the specification of the 
binormal in accordance with these requirements leaves us to consider the balance of forces in 
the plane of the curve. This situation is illustrated in Figure 2. The quantities v and —qz are 
the projections of the translational velocity v and body force q onto the plane of the curve. A 
planar curve is specifiable in terms of a single function 0{s), the tangential or Whewell angle, 
which we measure with respect to what we denote as the x direction, perpendicular to the z 
direction. We also specify a parameter (/>, the fixed angle between x and the projected velocity 
V. These elements are shown schematically in Figure 2. 



Figure 2: A planar catenary X with orthonormal frame (t,fi). The translational velocity and 
body force have been projected onto the t-n plane to obtain v and —qz, respectively. The flow 
Ft is unaffected. The parameter 0 is the fixed angle between the projected velocity v and the 
X direction (not labeled, but indicated by dotted lines) perpendicular to z. The body X can 
be described by the local angle 9 between its tangent t and x. 

Let us rescale quantities using an arbitrary reference length scale L and a velocity scale 
F= ■ V + + (q/ux)"^ , so that s —)■ s/L, |v| —>■ |v|/F, T —)■ T/F, q —)■ g/(z/j_F), and 

m —>■ mF/(z/_i_L). These rescalings lead to the constraint + v ■ v + = 1. They will need 
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to be revisited in Section 4.2, when we discuss a limiting case with vanishing v, for which z/j_ 
is meaningless. 

The planar, dimensionless form of (4) is 

mT^dsOh + (T + V ■ t) t + V ■ fin = dg (crt) — gz , (8) 

where a new parameter, the drag anisotropy ratio Ry = v\\/v±., appears. Using the relations 

/ t ^ / cose sine W X \ 

y n y y — sin e cos 6* y y z y ’ ^ 


we obtain the projections of (8) onto the two osculating frame vectors t and n, 

dgcc = Ry{T + V ■ X cos 6^ + V ■ z sin e) + q sin 6*, (10) 

adge = mT'^dse + (— v ■ x sin 6^ + v ■ z cos 6*) + g cos 6*. (11) 


The projections v-x and v-z appearing in these equations are constants related to the magnitude 
|v| and direction 0 of the translational velocity. The stress a can be eliminated to obtain a 
single second order equation for 6^, 


dg 


(— V ■ X sin 0 + V ■ z cos e) + q cos 6* 

^e 


(T + V ■ X cos 6* + V ■ z sin 6*) + g sin 6*. (12) 


Multiplying by the inverse of the bracketed quantity on the left, we obtain 


(—V ■ X sin 6^ -|- V ■ z cos e) + q cos e 


Ry(T -|- V ■ X cos 6^ -|- V ■ z sin 0) -|- g sin e 

dse 


(—V ■ X sin 6^ -|- V ■ z cos e) + q cos 6* 


(13) 

Both sides of this equation are total s-derivatives. Before proceeding, let us introduce a second 
rescaling, whose purpose is to simplify the expressions that will result from integrating the 
equation. We define A = a/ (v ■ x)^ + (v • z + g)^ , U = v ■ x/A, lU = v • z/A, T = RyT/A, 
Q = q/A, which leads to the constraint + {W + = 1. We may now write the integrated 

equation (13) as 


In 


where 


gjn. 

dse 

m 

m 


k 


C-f\n\f{e)\-{k-l)\n\g{e)\-{l-Ry)UQe, 

1 - U - (lU + Q) tan I 
1 + U + (lU + Q) tan I ’ 

—U sin e + (w + q) cos e , 

1 + RytJ^ + [RyW + Q){W + Q) , 


(14) 

(15) 

(16) 


and a term proportional to In A has been absorbed into the constant of integration C. Rear¬ 
ranging, we obtain a single first order equation for 0, 

9,0 = C'e(i-WO(^|/(0)f|g(0)|" . (17) 


Once e is determined by quadrature, a follows directly from (11) as 


a = mT'^ + 


Agje) 

dse 


(18) 


and the Cartesian position vector X follows by integration of the defining relation between 
tangent vector and tangential angle, = xcos0 -|- zsin6'. 
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The shape equation (17) is a hve-parameter family of dynamical systems, written in terms 
of six constrained parameters. A sixth independent parameter, the mass coefficient m, appears 
only in the stress equation (18) and does not affect the geometry of the solutions. We note, 
however, that the density of a massive cable acted on by gravity will indeed affect its geometry, 
through the body force parameter Q. Specifying the horizontal and vertical velocity parameters 
U and W is akin to specifying the translational velocity magnitude |v| and direction 0, and in 
the remainder of this paper we will use whichever description is simplest. The other quantities 
appearing are the axial flow parameter T, the drag anisotropy R,y, and a scaling constant C 
which is inherited from the limiting case of the classical catenaries. Due to the symmetry of 
the problem, we need only consider positive values of IJ and Q. We will also restrict Ry to the 
physically realistic range between zero (purely normal drag) and unity (isotropic drag). 

The right hand side of equation (17) is the product of an aperiodic exponential, a 27r- 
periodic function \f{9)\, and a 7r-periodic function \g{9)\. The 7r-periodicity of \g{9)\ is punc¬ 
tuated by zeroes at values 9o for which U tan 9o = W + Q. Excepting special cases, which 
will be discussed in Section 3.2, individual orbits generated by (17) will encounter either hxed 
points or poles at these locations, and cannot subtend an angle greater than tt. While it is not 
obvious by inspection, f{9) has as its distinguished points alternating zeroes and poles that 
coincide with the zeroes of g{9). A generic phase portrait, drawn in the 9-ds9 plane by varying 
the constant of integration C, will consist of a series of lobes separated by fixed points and/or 
poles. It turns out that a single lobe, subtending an angle of tt, suffices to describe all the 
possible generic shape equilibria. This is because a combination of lobe shift, sign change in 
axial velocity, and rescaling brings orbits in phase space into coincidence, 

_ ^ 

ds9 {9 o- 71 + 9, -f) = (r^) ’ (19) 

where 9q = arctan E [—f,f] , (20) 

and reflection about the 9 axis in phase space, achieved through a sign change in the scaling 
coefficient C —)■ —C, serves to identify configurations X in physical space. Details may be 
found in Appendix A. 

While the full form of equation (17) appears complicated, there are many situations in 
which it can simplify. For example, if U vanishes, so that the body’s translational velocity 
is collinear with the body force, the exponential term reduces to unity, and g{9) simplihes 
to ± cos 9. Similarly, if the axial flow T vanishes, so that the body translates rigidly, the 
complicated dependencies of f{9) become irrelevant. Curves for which both of these conditions 
hold can be found using an alternate method of solution. This is presented in Appendix B 
along with representations of the classical catenaries {k = 2) in 9-ds9 and K-dgK spaces. 

3 Solutions 

Before organizing our Endings in Section 3.2, we present several examples in Section 3.1 that 
represent the spectrum of typical solutions. 

3.1 Examples of generic shapes 

Almost all phase portraits generated by equation (17) are of three types, with either fixed 
points, poles, or alternating fixed points and poles occurring at intervals of tt. A fixed point 
solution is a straight line, with fixed angle 9o. All other orbits are qualitatively the same for 
each type of phase portrait. They are, respectively, heteroclinic connections, double-sided and 
single-sided blowup solutions, where the quantity that blows up is the curvature ds9 at finite 
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angle 6 and arc length s. Figure 3 shows four lobes from each type of phase portrait. The effect 
of the term exponential in 9 is apparent. 





heteroclinic connections 


double-sided blowup 


single-sided blowup 


Figure 3: The three generic types of phase portraits generated by equation (17), using C = 
±(1,2,3) and the following other parameters. Heteroclinic connections: 0 = —75°, |v| = 0.8, 
T = 0.5, q = 0.33, Ry = 0.5. Double-sided blowup: 0 = —88°, |v| = 0.8, T = 0.3, q = 0.52, 
Ry = 0.1. Single-sided blowup: 0 = —88°, |v| = 0.76, T = 0.4, q = 0.51, Ry = 0.4. See Figure 
2 and the adjacent text for explanations of the parameters. 


In this section, we present several examples of the three qualitatively different types of 
solutions. When displaying and comparing curves generated by different parameters, the scaling 
constant C is not particularly interesting. However, simply keeping C hxed does not always 
result in illustrative comparisons. We wish to compare several curves of equal length, while 
assuring display of a result corresponding to as much of the orbit as is possible. We achieve 
this by scaling solution curves such that a length of unity is associated with subtension of 
an angle slightly less than tt. This requires C to vary along with the other parameters. The 
condition we impose leads to a pair of curves corresponding to positive and negative values of 
C; one of these curves will be purely in tension (a > 0), and this is the curve we choose. Its 
partner, which may be partly or entirely in compression, and thus unstable, has a conhguration 
obtained by rotation of the tensile shape by tt radians. This representation is an efficient means 
of displaying families of curves, but often does not correspond to the variation of conditions 
arising naturally in a realistic boundary value problem. If this is not kept in mind, some of the 
curve comparisons may seem counterintuitive. 

Because orbits approach fixed points asymptotically at inhnite s, we cannot integrate the 
full TT radians and obtain hnite curves for display. Likewise, integrating up to a curvature 
singularity is numerically impractical. Instead, we choose a small cutoff angle e so that a 
typical shape will be represented as a curve of length unity subtending an angle of tt — 2e from 
Oq — + e to 6 q — e. The choice of cutoff affects the magnitudes of the curvature and stress, so 

care must be taken in comparing solutions with different values of e. All of the curves in this 
section (3.1) were generated with a cutoff e = 0.05. The normalization is achieved by using the 
fact that the length of the curve is 

( 21 ) 

1 C 
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so that a length of unity is achieved by setting 

C= / de^ (23) 

hrst substituting the expression obtained from (17) for the integrand on the right hand side. 
With this value of C, we integrate the system described in Section 2.2 from s = 0 to s = 1 or 
some other desired length. 

Each curve X is shown alongside its phase space orbit and its stress distribution. The 
orbits in the hgures will in general overlap, as they are extracted from different phase portraits. 
We set the mass coefficient m = 1 for all solutions in this paper, as it has no effect on the 
shapes, and simply shifts the stresses by a constant mT^. A dotted line on the end of a curve 
indicates that the curve could be continued out to inhnite length, while an open circle indicates 
that the curve terminates. Of course, one can always cut a solution subtending an angle less 
than TT from any of these mother curves. 

Figures 4-8 show catenaries rigidly moving at three different speeds in five different di¬ 
rections, with drag anisotropy taking the Stokes flow value = 0.5. These solutions are all 
heteroclinic connections corresponding to inhnitely long conhgurations. Not surprisingly, the 
alignment of the curves is set by competition between the body force —qz and translational 
velocity |v|v. There are also subtle differences in the conhgurations generated by a uniform 
body force or anisotropic drag forces, best seen in Figure 4 where the two vectors are aligned. 
The orbits in Figure 7 are signihcantly asymmetric, as is quite apparent from inspection of 
the conhgurations generated with a symmetric cutoh e from an inhnitely long mother curve. 
When looking at these hgures, it is important to recall that the solutions are generated by 
constraining the curves to be of unit length while subtending an angle just shy of tt radians. 
This is not equivalent to hxing the location of the curve ends and varying the velocity. The 
latter more realistic boundary value problem will be addressed in Section 4. The tension in all 
of these curves increases indehnitely as the curves go oh to inhnity in either direction. It is 
intuitively plausible that holding or towing an inhnitely long curve should require inhnite force, 
but in Section 3.2 we will show that this is not always the case. 

In Figure 8, the curves hip from being concave up to concave down. When v ■ z = W = 
—q, our scaling is pathological and the solutions are straight lines perpendicular to v and z. 
However, this is not the full story. The transition to concave down curves is not smooth. For 
a small window of parameters such that the velocity and body force are nearly balanced, we 
encounter double-sided blowup solutions corresponding to hnite-length curves, as in Figure 9. 
The tension prohle now peaks in the middle and, for these rigidly moving curves, goes to zero 
at the ends. Thus, these parachute-like objects are solutions to a free sedimentation problem 
for perfectly hexible hlaments. 

Figure 10 shows the ehect of a small amount of axial how on an inhnite length curve, with 
the ratio of drag to body forces maintained constant. Arrows indicate the direction of how. 
The biggest contribution to the variation in tension with axial how T is a centripetal shift of 
mT^. Figure 11 shows a hnite length curve with weak axial how, to be compared to the darkest 
curve in Figure 9 (while noting diherent overall scales). Here one should imagine two reels 
traveling downward and exchanging a length of cable between them. The tension at the ends 
is now mT^ rather than zero. 

With sufficiently large axial how, the solutions become of single-sided blowup type, cor¬ 
responding to semi-inhnite length curves. Figure 12 shows an example, with drag anisotropy 
Riy = 0.6. Here the tension is mT^ on the hnite end of the body, and goes to inhnity on the 
inhnite side. As will be discussed in Section 3.2, this is not always how the tension behaves in 
such curves. 
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Figure 4: Configurations (left), orbits (upper right), and stresses (lower right) for catenaries 
rigidly translating with velocity v under the action of a body force —gz. Darker curves indicate 
higher velocities. Parameters are 0 = 90°, |v| = (■ .1, 0.5, 0.95), T = 0, g = ( 0.87, 0.31), 

Ri, = 0.5, m = 1, e = 0.05. Infinite length curves generated by heteroclinic orbits. 
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Figure 6: As in prior figure, with 0 = 0°. 





Figure 7: As in prior figure, with 0 = —45°. 
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Figure 8: As in prior figure, with 0 = —90°. 





Figure 9: Finite length “parachute” catenaries. As in prior hgures, with parameters (j) = —90°, 
|v| = ( , 14, 0.75, 0.76), T = 0, g = ( .67, 0.66, 0.64), = 0.5, m = 1, e = 0.05. 
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Figure 10: Effect of axial flow T at fixed velocity to body force ratio \y\/q = 0.3. As in prior 
figures, with parameters 0 = —45°, |v| = (0.28, 0.29, 0.28), T =(-0.2, 0, 0.2), q = (0.94, 0.96, 
0.94), Ri, = 0.5, m = 1, e = 0.05. Infinite length curves. 





Figure 11: An axially flowing finite length catenary. As in prior figures, with parameters 
0 = -90°, |v| = 0.76, T = 0.13, q = 0.64, = 0.5, m = 1, e = 0.05. 
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Figure 12: An axially flowing semi-infinite length catenary. As in prior figures, with parameters 
0 = -90°, |v| = 0.78, T = 0.25, q = 0.57, = 0.6, m = 1, e = 0.05. 


3.2 Organization of solutions 

The transitions between these qualitatively different conhgurations are controlled by the signs 
and relative magnitudes of the exponents of f{9) and g{9) in equation (17). These exponents 
are the twice-rescaled axial flow T and the quantity k dehned in (16). Inequalities that govern 
the shapes are 


k + T > 0 and k — T > 0 

heteroclinic connections, 

(24) 

k + T < 0 and k — T < 0 

double-sided blowup, 

(25) 

else 

single-sided blowup. 

(26) 


The corresponding equalities are special cases that will be discussed below. 

There are also transitions in the behavior of the tension at the zeroes or poles. The extra 
factor of g{9) in equation (18) means that these do not coincide with the shape transitions. 
Instead, we have 


1-k + T >0 

and 

1-A;-T>0 

a —)■ mT^ 

on both sides, 

(27) 

1-k + f <0 

and 

1-k-f <0 

a —)■ cx) 

on both sides, 

(28) 



else 

a —)■ mT"^ 

on trailing side and 

(29) 




a —)■ cx) 

on leading side. 



Again, equalities are special cases that will be discussed below. 

Figure 13 is a simple planar map, spanned by the two exponents, that shows where different 
types of solutions are to be found. Teal, green, and orange regions denote heteroclinic, single¬ 
sided, and double-sided shapes, respectively. The purple and maroon borders are special cases 
for the shapes. The long- and short-dashed lines are special cases for the stress. Upon approach 
to a zero or pole, the stress approaches inhnity in the rightmost region of the diagram and a 
hnite value in the leftmost region. Note that there is a green strip of semi-inhnite curves with 


14 







finite tensions, and a teal strip of infinite curves with tension approaching a hnite value on 
one or both sides. These results seem strange. We would naively expect the dashed lines to 
coincide with the colored borders, so that the tension behavior echoed the shape of the curves. 

A more physically intuitive map may be drawn using the once-rescaled quantities. Solu¬ 
tions for shapes and stresses depend on six independent parameters, but the scaling parameter 
C and mass coefficient m do not affect our qualitative classihcation. Fixing the angle of trans¬ 
lation 0 and the drag anisotropy Ru leaves us with two independent parameters. These may be 
represented as coordinates on a 2-sphere in the space spanned by the magnitudes of the axial 
flow, translational velocity, and body force, T, |v|, and q, related by -|- v ■ v -|-= 1. Due to 
the symmetry of the problem, we need only consider one octant of the sphere. Figure 14 shows 
an example for 0 = —88° and = 0.5, with colors and dashing retaining the same meaning 
as in Figure 13. Solutions on the poleward T side of the dashed lines have tension inhnite 
on their leading ends and hnite on their trailing ends. Those in the teal heteroclinic regions 
that include the |v| and q corners have inhnite tension on both ends, while those in the region 
that entirely encloses the orange double-sided blowup region have hnite tension on both ends. 
Moving along a lattitude line keeps the axial how hxed while varying the relative strengths of 
the translational drag forces and body force. Moving along a meridian keeps the ratio of drag to 
body forces hxed while varying the axial how. For ease of illustration and viewing, we use the 
orthographic projection of the octant onto the |v|-g plane, a circular sector shown in Figure 15. 
Physical chemists or the metallurgically inclined may note that this is somewhat reminiscent of 
a ternary phase diagram. The corners of the sector represent the pure classical catenary (g), a 
neutrally buoyant towed cable (|v|), and a lariat snaking along itself (T). However, the latter 
limit, corresponding to the pole of the sphere, has no meaning here because the rescaling by A 
is equivalent to dividing by zero. Similarly, when the angle of translation is —90°, the rescaling 
breaks down on the central meridian where g and |v| are equal. 

f 


k 

Figure 13: A planar map of solution types, spanned by two exponents k and T that appear 
in the shape equation (17). The teal, green, and orange regions correspond to heteroclinic 
connections, single-sided blow up, and double-sided blow up solutions, respectively. The purple 
border between teal and green, and the maroon border between orange and green, correspond 
to special cases of the shape equation. The long- and short-dashed lines correspond to special 
cases of the stress equation (18). With respect to the dashed lines, the asymptotic end stresses 
are inhnity in the region on the right, mT^ in the region on the left, and one of each in the 
region in between. 
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Figure 14: A spherical map of solution types spanned by the magnitudes of the translational 
velocity |v|, axial flow T, and body force q. Fixed parameters are 0 = —88° and = 0.5. As 
in Figure 13, colored regions correspond to heteroclinic connections, single-sided blow up, and 
double-sided blow up solutions. With respect to the dashed lines, the asymptotic end stresses 
are inhnity in the regions adjacent to the |v| and q corners, mT^ in the region enclosed by the 
long-dashed line and the equator, and one of each in the region adjacent to the T pole. 



Figure 15: Orthographic projection of the octant of the (j) = —88°, Ry = 0.5 sphere from Figure 
14 onto the |v|-g plane. 
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In Figure 16, we show five curves taken from the T = 0.2 latitude line on the sphere for 
nearly downward translation 0 = —89.9° and Stokes drag anisotropy Ry = 0.5, and correspond¬ 
ing orbits in a dvr-wide slice of phase space. The teal, green, and orange curves correspond to 
the heteroclinic, single-sided, and double-sided solutions already discussed. In between these 
are purple and maroon curves that subtend an angle of 27r rather than vr. They are special 
cases corresponding to the boundaries between the generic regions. When k = ±f, f{9) and 
g{9) combine to form 27r-periodic functions. New orbits and configurations are formed from 
the fusion of two heteroclinic or two double-sided blowup solutions. In the figure, we have 
scaled configurations as before, with a length of one unit per angle of vr subtended, but used 
the scaling C = 4(1 — Ry)UQ7r to draw the orbits. This scaling prevents overlap of the orbits 
in the range shown, but is insufficient to prevent overlap at more negative values of 9. The 
result is a set of curves that seem to fit together within the range, and are suggestive of a phase 
portrait. It is unknown whether the phase space, or some slice of it, can be foliated by orbits 
generated by all the intervening points along our circumferential cut through parameter space. 
If so, it suggests that there is an additional symmetry hidden in our system. 

Figures 17 and 18 show examples of 27r-subtending special cases for perfectly downward 
translation 0 = —90°. These two curves are symmetric about their midpoint, being formed 
from a vr-subtending curve and its mirror image. Their tensions are symmetric about the value 
mT^. Because the tension takes on negative values in portions of the concave-upward portions 
of the curves, these solutions are unstable. Note that Figure 17 was generated with Ry = 0.4, 
and that both curves were generated with a smaller value of the cutoff, e = 0.001, than in 
Section 3.1. These curves have been scaled to have a length of two units. 




Figure 16: Five axially flowing catenaries taken from the dotted locations along the T = 0.2 
latitude line in the projection (center) of a 0 = —89.9°, Ry = 0.5 spherical map octant. The 
configurations (right), and several periods of orbits (left), correspond to a generic vr-subtending 
heteroclinic connection, a 27r-subtending heteroclinic special case, a generic vr-subtending single¬ 
sided blow up solution, a 2vr-subtending double-sided blow up special case, and a generic vr- 
subtending double-sided blow up solution. Configurations are drawn by adjusting C so that a 
length of unity corresponds with a subtended angle of approximately vr, while orbits are drawn 
using the scaling C = 4(1 — Ry)UQ'n' to eliminate overlap in the range shown. The other 
parameters are |v| = (0.83, 0.81, 0.78, 0.76, 0.74), q = (0.52, 0.55, 0.59, 0.62, 0.64). 
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Figure 17: An axially flowing, self-intersecting loop de loop catenary. Details as in Figures 
4-12, with parameters 0 = —90°, |v| = 0.82, T = 0.2, q = 0.53, Ri, = 0.4, m = 1, e = 0.001. 
Infinite length curve subtending 27r due to the special condition k = T. 





e 


Figure 18: An axially flowing C-clamp catenary. As in prior figure, with parameters 0 = —90°, 
|v| = 0.78, T = 0.1, q = 0.61, = 0.5, m = 1, e = 0.001. Finite length curve subtending 27r 

due to the special condition k = —T. 
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Another curious special case is found at either of two triple points, where the three regions 
meet. This occurs on the equator, where T = 0 and the curve rigidly translates. In Figure 19, 
we show three curves taken from the T = 0 latitude line for a translational angle 0 = —88.5° 
and Stokes drag anisotropy Ry = 0.5. At the triple points, both exponents T and k are zero, 
and the solutions are logarithmic spirals ds9 = . The tension is a damped sinusoid, 

and the compressive regions make the curves unstable. A portion of spiral subtending 127r, 
scaled to be 12 units long, is shown in Figure 20. No cutoff e is necessary to generate this 
curve, as all the poles and hxed points have disappeared. For perfectly downward translation 
0 = —90°, one triple point spiral becomes a straight line with v = —qz, and the other becomes 
a circle dgO = C with v = — and an undamped sinusoidal tension. To our knowledge, 

the circle is the only periodic solution, excepting the arbitrary shapes (mentioned in Section 
2.1) that are possible for isotropic drag Ry = 1. It is also a freely sedimenting solution. It is 
perhaps worth pointing out that, unlike elastic curves with bending resistance, and with the 
exception of = 1 arbitrary shapes, all of our augmented catenaries are non-inflectional. 

Examples of stress transitions upon crossing dashed lines are shown in Figures 21 and 22. 
The solutions are generated with Ry = 0.6 and a small value of the cutoff e = 0.005. While 
the conhgurations and orbits display no qualitative changes, the tension on one side of the 
curve changes from asymptotically inhnite to the hnite value mT^. For the special cases on the 
boundaries, the tension takes on an intermediate hnite value. In these hgures, we encounter 
the unexpected phenomena of inhnitely long curves with hnite tension on one side, and semi- 
inhnite curves with hnite tension on both sides. There also exist cases, not shown, of inhnitely 
long curves with hnite tension on both sides. 

While this behavior is not fully understood, we may gather a clue to its origin by consid¬ 
ering rigid collinear translations (T = = 0) in the part of the heteroclinic region with hnite 

tension. As is evident from inspection of Figure 13, these solutions occur when 0 < /c < 1. 
Equation (16), the constraints on twice-rescaled quantities, and the positivity of Q then indi¬ 
cate that RyW -|-Q > 0 while W + Q = —1. In terms of the once-rescaled quantities, Ry\'v\ < q 
and |v| > q. This means that portions of the curve aligned with the direction of translation 
and body force are moving slower than they would be if allowed to free-fall at their terminal 
velocity, while portions perpendicular are moving faster than their terminal velocity. 

Finally, we show the ehects on the spherical map of varying the two remaining relevant 
parameters. Figure 23 shows maps for six angles of translation at hxed Stokes anisotropy 
Ry = 0.5. The double-sided blowup region is only present for sufficiently aligned translation. 
As the direction of translation approaches that of the body force, the dotted stress transition line 
collides with the equator at a point, which then splits into two spreading points on either side 
of a newly formed hnite-stress region. Likewise, the green single-sided blowup region collides 
at a diherent point, its contacts then spreading around a new double-sided blowup region. No 
changes occur on the meridians where |v| or g are zero, as the quantity 0 has no meaning 
there. Figure 24 shows maps for six diherent drag anisotropies at hxed angle of translation 
0 = —85°. The double-sided blowup and counterintuitive inhnite-length, hnite-stress region 
expand with increasing anisotropy; the latter occupies half of the sector when Ry = 0. At this 
limit, there is no tangential drag, the axial how parameter T has no ehect on the shape, and 
the entire sector is radially symmetric. Both regions disappear as Ry approaches unity and the 
anisotropy becomes sufficiently small. Note that for the isotropic Ry = 1 map in this hgure, 
the single-sided region and stress transition line do not touch the equator. The stress transition 
on the g = 0 meridian occurs when the axial how T and translational velocity magnitude |v| 
are equal, independently of the value of Ry. At certain limits, the single-sided blowup region 
becomes a perfect ellipse, the orthographic projection of a region between two small circles on 
the sphere. When 0 = —90°, the ellipse is given by [1-|- (l-|-i?i.)^]|vp-|-4(l-|-i?,^)|v|g-|-5g^ = 1, 
and when Ry = 1, it is given by 5|vp -|- 8|v|gsin0 -|- 5g^ = 1. 
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Figure 19: Three rigidly translating catenaries taken from the dotted locations along the T = 0 
equator in the projection (center) of a 0 = —88.5°, = 0.5 spherical map octant. The 

conhgurations (right), and several periods of orbits (left), correspond to a generic vr-subtending 
heteroclinic connection, a logarithmic spiral special case, and a generic vr-subtending double¬ 
sided blow up solution. Configurations are drawn by adjusting C so that a length of unity 
corresponds with a subtended angle of approximately vr, while orbits are drawn using the 
scaling C = 4(1 — R^)IJQti to eliminate overlap in the range shown. The other parameters are 
|v| = (0.82, 0.8, 0.78), q = (0.57, 0.6, 0.63). 





9 

Figure 20: A logarithmic spiral with damped sinusoidal stress is a rigid catenary. As in Figures 
4-12 and 17-18, with parameters (p = —88.5°, |v| = 0.79, T = 0, g = 0.61, R^ = 0.5, m = 1, 
e = 0. Rigidly translating infinite length curve subtending infinite angle due to the special 
condition k = T = 0. 
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Figure 21: Projection and configurations (left), orbits (upper right), and stresses (lower right) 
for three single-sided blowup catenaries taken from the dotted locations along the T = 0.25 
latitude line, with hxed parameters 0 = —89° and R^, = 0.6. Darker curves indicate higher 
velocities. In crossing the long-dashed boundary line, the stress on the leading side of the body 
transitions from infinite to the hm value mT^, with the special case on the boundary taking 
on an intermediate hnite value. The other parameters are |v| = ( , 0.76, 0.78), q = { , 

0.60, 0.57), m = 1, e = 0.005. 
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Figure 22: As in prior hgure, for three heteroclinic catenaries taken from the dotted locations 
along the T = 0.2 latitude line, with |v| = ( 0.88, 0.93), q = ( 0.41, 0.30). In 

crossing the short-dashed boundary line, the stress on the trailing side of the body transitions 
from infinite to the hm value mT^, with the special case on the boundary taking on an 
intermediate hnite value. 
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(a) (t) = 90° 


(b) (/- = 0° 


(c) 0 = -60° 





(d) (/. = -81.79° (e) (/) = -88° (f) </> = -90° 

Figure 23: Effect of translational angle 0 on spherical maps with hxed drag anisotropy Ru = 0.5, 
shown through orthographic projections onto the |v| — g plane. Colored regions, colored lines, 
and dashed lines have the same meanings as in Figure 14. 
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(a) i?i/=0 


(b) Ru=0.l 


(c) R^=0.2 



(d) i?^=0.4 


(e) R^=0.6 


(f) Ru=l 


Figure 24: Effect of drag anisotropy Ry on spherical maps with hxed translational angle 0 = 
—85°, shown through orthographic projections onto the |v| — g plane. Colored regions, colored 
lines, and dashed lines have the same meanings as in Figure 14. The solid and dashed lines in 
the bottom right panel do not touch the T = 0 equator. 
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4 Boundary value problems 

So far, we have attempted to compactly display a multiparameter family of curves using bulk 
physical quantities such as velocities, and geometric variables such as angles. These quantities 
are only implicitly related to practical boundary conditions, which usnally depend on positions, 
or to global constraints such as the total length of a cable. In this section, we discuss a 
few boundary valne problems of potential interest. The hrst example involves a fnlly three- 
dimensional solution, where we recognize that our planar curves and the forces acting on them 
are embedded in a three-dimensional world. The last example lives in the limit of vanishing 
translational velocity, and so reqnires a new choice of rescaling. 

As described below in Section 4.1, a three-dimensional problem is first rednced to a two- 
dimensional problem corresponding to the system described in Section 2.2. Two-dimensional 
problems are solved using two-parameter Newton shooting. In the following examples, the 
curves have specihed endpoint positions and parameters, with the exception of the scaling 
parameter C. The system described in Section 2.2 is integrated over the appropriate arc length 
s from one end with a pair of given Cartesian coordinates, shooting for the coordinate pair of 
the other end by adjusting C and the initial angle 6. 

4.1 Towing a cable from two points 

We retnrn to dimensionful quantities. Consider a fixed length of cable, acted on by a body 
force q, being towed with velocity v by two spools separated by a vector R. In general, these 
are non-orthogonal vectors spanning three dimensions, and onr hrst task is to determine the 
plane containing the cable. This plane is independent of any axial how that might be present. 
We will identify it by identifying the nniform binormal vector b of the curve. For this we use 
the force balance on the binormals, eqnation (7), and the two additional equations b ■ b = 1 
and b ■ R = 0. Once the binormal is known, we may dehne the projections of the translational 
velocity and body force as v = v — v ■ bb and —gz = q — q ■ bb. The translational angle 0 is 
snch that v ■ x = |v| cos 0, with x specihed by the relations x ■ b = 0, x ■ z = 0, and x ■ x = 1. 
Nondimensionalizing leads to the qnantities appearing in eqnation (8). 

Three cnrves, representing three states of a single cable, are shown in Figure 25. The 
curves are rigidly towed in the same direction at diherent speeds, are identical in length, and 
their ends are separated by the same vector R. For this example, R and the velocity v are 
perpendicular to the body force q, but not to each other. As the towing velocity increases, 
the curves pivot to become more aligned with the velocity, and also change shape within their 
planes. 

The ehect of drag anisotropy is shown in Figure 26 for the case of pure rigid towing in the 
absence of body force. The problem is inherently two-dimensional for these nentrally buoyant 
cables. It is also eqnivalent to the collinear case discussed in Appendix B. The shape equation 
takes the particularly simple form dgO = C* cos*'^’*'^''^ 6^ for positive cos6'. For isotropic drag 
Ry = 1 the solutions are eqnivalent to the classical catenaries; the translational velocity has 
the same effect as a body force pointing in the opposite direction. Another case amenable to 
farther integration is that of pnrely normal drag = 0, where the arc length may be implicitly 
dehned throngh s = A In . 
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Figure 25: Two views of oblique towing of a unit length cable in a direction perpendicular 
to the body force. Darker curves indicate higher velocities. Unit vectors in the directions of 
the body force q, translational velocity v, and end-to-end separation R are shown. The angle 
between v and R is 45°, and they are both perpendicular to q. Dimensionful parameters are 
V = (', 8, 12), |q| = 4, z/|| = 0.4, u± = 0.8, and |R| = 0.8. 



Figure 26: Oblique towing of neutrally buoyant cables with different drag anisotropies Ry =( , 
0.5, 1). Darker curves indicate lower anisotropy. The dimensionless length of the cables is 
unity, and the end-to-end separation magnitude |R| = 0.8. 


4.2 Reeling a cable between two points 

Consider an axially flowing catenary reeled between hxed supports. In this limit, is mean¬ 
ingless, so we must rescale the quantities in equation (4) differently. Let us use an arbitrary 

reference length scale L and a velocity scale F = ; so that s —)■ s/L, |v| —)■ |v|/F, 

T — )■ T/F, q —)■ g/(z/|jF), andm —)■ mF/(z/||L). These rescalings lead to the constraint T^-|-g^ = 1 
on the axial flow and body force. Recall from Section 2.1 that if the body force q vanishes, 
then the only equilibrium solution is a straight line. 

Following the same steps as in Section 2 leads to the shape and stress equations 


d.9 = C 


tan 


e 


1 -|- tan I 


a = mT^ + 


e 
2 

q cos 9 


cos^ 9 , 


(30) 

(31) 
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We note in passing that essentially the same derivation can be performed if the drag is any 
nonlinear function of the axial velocity. This is relevant to the “lariat chain” conhguration of 
a loop of string hanging from a rotating wheel [25, 30, 31]. 

Due to the constraint, the dynamical system depends on only two parameters, the scaling 
constant C and the ratio T/q indicating the relative importance of the axial flow. The transition 
from inhnite to semi-inhnite stress occurs at T/q = 1, while that from heteroclinic to single¬ 
sided phase portraits occurs at T/q = 2. 

Figure 27 shows the effect of increasing axial flow on the curve conhgurations. The ca¬ 
ble lifts and, somewhat counterintuitively, shifts towards its trailing end. These qualitative 
effects can be observed in real “lariat chains”. Figure 28 shows two curves of different lengths 
corresponding to the special 27r-subtending heteroclinic case T/g = 2 at the purple boundary 
between heteroclinic and single-sided phase portraits. Unlike the generic solutions in Figure 27, 
the curves in Figure 28 can subtend an angle greater than tt. Note that a portion of the longer 
cable is in compression, so this solution is unstable. However, increasing the mass coefficient 
m could bring the entire cable into tension. 



Figure 27: Reeling of a unit length cable at different speeds, with end-to-end separation magni¬ 
tude |R| = 0.92. Darker curves indicate faster reeling. Parameters, rescaled using the tangential 
drag coefficient i/y, are T = ( , 0.88, 0.931) and q = { , 0.47, 0.365). 




Figure 28: Conhgurations (left) and stresses (right) for curves cut from the 27r-subtending 
heteroclinic special cases of the reeling problem. The ratio of axial how speed to body force is 
T/q = 2, the ratios of length L to end-to-end separation magnitude are T/|R| = (2, 4), and 
the mass coefficient m = 1. 


5 Discussion 

Our study has brushed against several established topics in the mechanics of huid-hexible 
structure interaction, and raised a few new questions. We briehy discuss these issues here. 

Our analysis allows us to generate a collection of spherical maps, exemplihed by Figure 15, 
and begin to think about what paths through this terrain might correspond to relevant boundary 
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value problems or interesting bifurcations. Bifurcations that create or destroy solutions likely 
tell us something about stability, a topic we have not explored in this work. We did not 
encounter such phenomena in our incomplete survey of boundary value problems in Section 4, 
but they might be found somewhere in the hve-dimensional parameter space of conhgurations. 
The six-dimensional space of solutions includes bifurcations in the asymptotic behavior of the 
tension. It remains to be seen whether such transitions have any signihcant effect on realistic 
physical situations when the cables are of hnite, indeed moderate, length. 

With respect to paths through parameter space. Figure 16 is particularly intriguing and 
raises an interesting question. Given a family of dynamical systems that depend on several 
parameters, when can one generate a phase portrait by following a path through parameter 
space? The paths we found to generate Figures 16 and 19 do not change the location of hxed 
points and poles, but do involve bifurcations that convert one type into the other. The ends 
of the path he inside the parameter space rather than on one of its boundaries. We have not 
established that these paths generate true phase portraits. If they do, it implies that there exists 
a physically or geometrically relevant parameter not captured by our analysis, and implicitly 
related to the parameters we chose for representation of our system. 

The free sedimentation of a hlament is a limiting case of our system. This problem’s 
history begins with the question of motion and reorientation of DNA or other macromolecules 
under centrifugation or an electric held [32-35]. Many have studied the sedimentation of rel¬ 
atively stiff hbers [36-39] and hber suspensions [28, 40] at low Reynolds number. Unlike the 
generic catenaries, these shapes can bend through angles greater than vr [38, 39]. One simula¬ 
tion demonstrated that at sufficiently low bending stiffness, the system bifurcates from steady 
sedimentation to periodic motion [38]. In the limit of highly hexible thermal chains [41, 42], 
the behavior appears to be chaotic, resembling the recirculating blob and escaping tail formed 
by clouds of falling particles at low to moderate Reynolds numbers [43, 44]. Our hndings are 
consistent with chaotic behavior for an athermal system in the perfectly hexible limit, in that 
there exists a multiparameter continuum of shapes and terminal velocities corresponding to the 
double-sided blowup solutions of Section 3.2. We would expect transient or terminal chaos if 
some or all of these are unstable. Note, however, that any bending stihness would not only reg¬ 
ularize the curvature singularities at the ends of the hlament, but require vanishing curvature 
there to satisfy a moment balance boundary condition. Hence, inferences about even weakly 
stih hlaments from results in our singular limit may not be well advised. On the other hand, 
for an elastic rod— as opposed to a macromolecule with elastic properties that could be more 
exotic than those of a simple rod— at low Reynolds number, nonlocal hydrodynamic forces 
will dominate over bending. These next-order interactions are likely responsible for much of 
the nonlinear phenomena observed in the literature. They are also likely to destabilize many, if 
not all, of the sedimentation solutions found in the present study. All this aside, the nature of 
the transition from steady to periodic to presumably aperiodic sedimentation with decreasing 
bending stihness is a topic worthy of further investigation. 

Extending the towing results of Section 4.1 to steady rotations and the signed quadratic 
drag forces rehective of high Reynolds number would aid in understanding towed cable and 
drogue dynamics in the air or ocean. In the context of single-end towing problems, much 
attention has been paid to linear and circular maneuvers, as well as combinations of these basic 
elements to form hairpin turns [45-67]. Equilibrium shapes and tensions are of interest, as 
are small vibrations superimposed on these background solutions. Pode captured the salient 
features of the planar rigid towing problem in terms of four quadratures [45]. Building on 
this work, Eames [46] provided analytical expressions for the Cartesian position vector of a 
cable towing a heavy drogue and translating perpendicularly to the body force, for several 
diherent drag assumptions including isotropic linear drag. He also discusses other interesting 
boundary value problems. While the non-analyticity of signed quadratic drag is inconvenient. 
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it is a surmountable difficulty, and it appears upon preliminary glance that these equilibria 
can be reduced to either first- or second-order dynamical systems. Such representations would 
be of lower order than the systems of equations in some of the preceding references, and of 
considerable advantage, in the search for equilibria, over the full-blown numerical schemes in 
the others. 

In the limiting case of reeling discussed in Section 4.2, any nonlinear relationship between 
velocity and drag force can be accommodated in a shape equation much like (30). Human- 
scale experiments using chains or ropes in a “lariat chain” configuration [25], for which we 
expect the drag to be quadratic in form, should afford a test of the theory. Another theoretical 
study of this configuration may be found in an early paper by Gregory [68], whose results are 
consistent with our shape bifurcation from heteroclinic to single-end blowup when drag forces 
become sufficiently large. It should also be possible to experimentally search for configurations 
corresponding to transitions between solution types. At a particular speed, the cable should be 
able to subtend an angle greater than vr, though for a given length the stability, and therefore 
the observability, of this solution will depend on the weight of the material. The existence of 
such a solution only under such restrictive conditions on speed is surprising, and this result will 
remain suspect until it is observed experimentally. 

The problem of laying and recovery of cable is a hybrid of towing and reeling problems, 
with a history that begins with the work of Thomson [6] and Airy [5] in the mid-19th century, 
and continues into the modern day [69-72]. Thomson neglected drag forces, while Airy unreal¬ 
istically assumed the forces to be linear and isotropic, corresponding to our limit = 1. The 
mid-20th century work of Zajac [69] is particularly worthy of note in the light of the present 
study in that, with particular assumptions on the configurations and boundary conditions, he 
considered the two-dimensional translation problem with axial flow and quadratic drag as a 
pair of coupled first order ordinary differential equations akin to our equations (10) and (11). 

Extension of this work to non-uniform high Reynolds number flows, and the non-planar 
shapes they can induce, is relevant to sonar, cable laying, mooring lines, and newer applications 
such as tethered wind energy devices [71-74]. It may also have some biological relevance. The 
deformation of plant life in response to non-uniform winds and waves, euphemized as “reconfig¬ 
uration” by Vogel [75], is a topic attracting much recent activity in the form of experimental and 
modeling studies of living and inanimate blade-like planar structures [12, 13, 76-83]. Possible 
complications here include bending resistance, and non-uniform buoyancy due to variation in 
morphology or health along the blade [77, 80]. In some of these, as well as other, high-Reynolds 
number systems, it may be reasonable and analytically helpful to neglect tangential drag forces. 

From a more fundamental point of view, the question of rotating, axially flowing equilibria 
subjected to linear drag forces would provide a simple framework to study the interaction 
between gyroscopic and dissipative mechanisms for breaking the axial flow symmetry of a one¬ 
dimensional continuum living in three dimensions. 

6 Summary 

We have presented and organized a six-parameter family of solutions, containing a five-parameter 
family of shapes, that extend the classical problem of the catenary to include linear drag forces 
resulting from translation and axial motion. The shapes are described by a first-order dynam¬ 
ical system. Limiting cases include sedimentation, towing, and reeling of cables. The effects of 
variation of parameters on solutions were displayed, first in a comprehensive but abstract way, 
and then with reference to a few specific problems with physically realistic boundary conditions 
and an integral constraint on length. These solutions may serve as a starting point for more 
advanced studies of flexible structure equilibria in fluid flows. 
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Appendix A Symmetries of the shape equation 


We wish to show that the two functions 


fiO) 


1 - f/ - (W + Q) tan I 
1 + f/ + (W + g) tan I 
—U sin 9 + {W + g) cos 9 


have either a pole or a zero at the same values of 9. The function g{9) has zeroes at angles 
arctan , spaced vr apart. The denominator or numerator in f{9) will vanish at angles 
2arctan We use the identities arctan a: = 2 arctan and arctan a: + arctan - = 

and the constraint U'^ + (W + g) =1, which can also be written as to conclude 

that arctan = 2 arctan = 2 arctan = vr — 2 arctan ■ 

We also wish to justify the rescaling of 


ds9 = \f{9)f \g{9)\^ 


after a shift of vr and change in sign of T, 

d,9 {9o-n + 9, -f) = (l^) > 

where tan | = ^ = ^ . 

The function \g{9)\ is vr-periodic, the behavior of the exponential under the transformation is 
simple, and the only place where T appears is in the exponent of the function \f{9)\. We need 
only demonstrate what happens to this function. We use the identities tan (x + y) = 

and tan (x ± |) = —cot x to hnd that f{9o + 9) = tan | and f{9o — tt + 9) = cot |, 

and therefore 

[/(e„-,r + «)]-^= ([^]) !/(«„ +«)f. (32) 

Finally, the constraint on U, W, and Q implies that f/ < 1, so |f/ — l| = 1 — U. 

The transformations and scaling are displayed graphically in Figure 29 for an exemplary 
choice of parameters. 
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Figure 29: The effect of T on the function dgO. Two roots are labeled with black dots. The 
parameters are 0 = —85°, |v| = 0.5, T =(-0.4, 0, 0.4), q = (0.77, 0.87, 0.77), C = 3. The 
transformation (19) and additional prefactor scaling relates a blue orbit to a red orbit in an 
adjacent lobe. The corresponding reflected red orbit (not shown) below the 9 axis generates 
the same physical conhguration as the original blue orbit. 


Appendix B Alternate derivation and representation of 

collinear systems 


Here we consider rigid translations collinear with the body force, including as limiting cases the 
classical catenaries and neutrally buoyant towed cables. This section makes use of the relations 
(1) that describe derivatives of the frame vectors. 


For rigid translations, T = 0, and the projections (5-7) 
problem onto the frame vectors are modihed accordingly. 

of the dimensionful three-dimensional 

z/||V ■ t = dsa -h q ■ t , 

(33) 

i/_i_v ■ h = (TK -|- q ■ n , 

(34) 

iy±-v ■ b = q ■ b . 

(35) 


We may note immediately that combining (34) and an s-derivative of (35) leads to the equation 
aKT = 0, implying that any solutions with nontrivial stress a and curvature k must be planar 
(r = 0). Therefore, we restrict consideration to the planar problem obtained by projecting the 
rigid form of the dimensionless equation (8) onto t and n. We also let v = nz, collinear with 
the body force, so that 

{RyV + q)i -t = dgU , (36) 

{v + q)z ■ Yi = an . (37) 

Take a derivative of (36) and use (37), and take a derivative of (37) and use (36), to obtain 
two equations, 

dla — = 0 , 
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(38) 

(39) 





depending on a single parameter Ay = ■ This parameter is equal to Ry for neutrally 

buoyant towed cables (g = 0), and equal to unity for the classical catenaries (n = 0) and for 
the equivalent case of isotropic drag {Ry = 1). Divide (39) by ctk to get, for positive a and n, 

9s In (T + 9s In K = 0 , 

na = B , (40) 

B constant, then insert (40) into (38) to get 

_ / 2+Aty \ 

9sV - AyB'^a V y = 0 , (41) 

and hnally multiply by ds<J and integrate to get 

{dsaf + AlB^a-^^ = Cy . (42) 

Cy constant. If the stress a is obtained from equation (42), the Cesaro equation of the curve 
follows as 

_ ( l+Ay \ 

K = Ba y . (43) 

For many values of Ay^ the analytical solution of equation (42) is implicit and complicated. 
However, when Ay = 1, corresponding to the classical catenaries or isotropic drag, the stress 
is a square root of a quadratic polynomial in s, and the curvature k is the reciprocal of that 
polynomial. This is a well known representation of the catenaries. For general values of Ay, 
the curve shape is more easily obtained by inserting (40) in (42) and multiplying through to 
eliminate denominators, 

2(2Ay + l) / ^ \ 

{dsK) + {Ay + 1) K Ay + 1 t^Ay + l - Cy] = 0 . (44) 


The dynamical system (44) can be derived from the rigid collinear form of (17), dgO = C cos^ 6 

2 

for positive cos^. We note that k = 1 + Ay and identify Cy = C^+^^. In Figure 30 we display 
two representations of the classical catenaries, for which Cy = C, as heteroclinic orbits in O-dgO 


space and as homoclinic orbits in K-dgK space. 



Figure 30: Two phase space representations 
curvature variable. On the left, ds9 = C cos' 
a = 0 = ±(1,2, 3). 



K 


of the classical catenaries, using an angle or 
6, on the right, (9*^)^ = 4A {Cy — k), using 
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